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ABSTRACT 

We study the influence of initial conditions on the magnetic field amplification during the 
collapse of a magnetised gas cloud. We focus on the dependence of the growth and satura- 
tion level of the dynamo generated field on the turbulent properties of the collapsing cloud. 
In particular, we explore the effect of varying the initial strength and injection scale of turbu- 
lence and the initial uniform rotation of the collapsing magnetised cloud. In order to follow 
the evolution of the magnetic field in both the kinematic and the nonlinear regime, we choose 
an initial field strength of 1 /iG with the magnetic to kinetic energy ratio, E^/E\^ ~ 10~^. 
Both gravitational compression and the small-scale dynamo initially amplify the magnetic 
field. Further into the evolution, the dynamo-generated magnetic field saturates but the total 
magnetic field continues to grow because of compression. The saturation of the small-scale 
dynamo is marked by a change in the slope of Bj p^l'^ and by a shift in the peak of the mag- 
netic energy spectrum from small scales to larger scales. For the range of initial Mach numbers 
explored in this study, the dynamo growth rate increases as the Mach number increases from 
^rms/cs ~ 0.2 to 0.4 and then starts decreasing from ^rms/cs ~ 1.0. We obtain saturation val- 
ues of E^/E\^ = 0.2 — 0.3 for these runs. Simulations with different initial injection scales 
of turbulence also show saturation at similar levels. For runs with different initial rotation of 
the cloud, the magnetic energy saturates at E^/E\^ ~ 0.2 — 0.4 of the equipartition value. 
The overall saturation level of the magnetic energy, obtained by varying the initial conditions 
are in agreement with previous analytical and numerical studies of small-scale dynamo action 
where turbulence is driven by an external forcing instead of gravitational collapse. 

Key words: stars: formation - methods: numerical - magnetic fields - turbulence. 



1 INTRODUCTION 

Magnetic fields are ubiquitous in astrophysical systems and their 
study forms an active area of research today. Radio observations 
over the last few decades have revealed that galaxies and galaxy 
clusters host magnetic fields. The total magnetic field strength in 
nearby spiral galaxies is 9 — 15 /xG { [Beck & Hoernes|[l996[ 
|Beck|2004j ) while observations of cluster magnetic fields show that 
fields are at the /xG level, with values up to tens of /iG at the 
center of cooling core clusters ( Carilli & Taylor 2002', 'Go voni &| 
|Feretti|[T004| ). Recent observations also point to the existence of 
magnetic fields in the high-redshift universe ( [Bernet et al.|2008| ). 
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One plausible mechanism of the origin of such magnetic fields is 
the dynamo process where energy in the turbulent fluid motions is 
tapped to amplify the magnetic field. Turbulence is ubiquitous in all 
astrophysical systems ranging from protostellar accretion disks in 
individual star forming clouds to the interstellar medium (ISM) in 
galaxies and possibly also in the gaseous media of galaxy clusters 
and cosmic filaments. In fluid dynamics, turbulence is described as 
a flow regime characterised by chaotic motions and involves the 
cascade of energy from the scale of the largest eddy to the smallest 
scale eddy. In particular, the small-scale dynamo proc ess ([Kazant- 



sev 1968 ; Vainshtein 1982 ; Schekochihin et al. 2002','Boldyrev & 
Cattaneo 2004, Sur et al. 2010, Federrath et al. 2011b, Schober 
et al.|20lT] ) can lead to rapid amplification of initial seed magnetic 



fields (see 



Brandenburg & Subramanian|2005| for a review). Ear- 
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lier work by 'Beck et al.| ( |199 7) proposed that such dynamo action 
is responsible for producing seed magnetic fields for the galactic 
large-scale dynamo. The same mechanism is found to amplify mag- 
netic fields in galaxy clusters ( ^Dolag et al.|1999|[Subramanian et al.| 
[20061 |Xu et aL||2059l[20rT]) and also in the cosmological large- 
scale structure ("R yu et al.||2008| ). Mergers during cluster forma- 
tion can also lead to turbulence and intense random vortical flows 
([Schindler & Mueller] 19 93 , KulsmdetaL|l997 , Bryan & Norman 
1998[ IMiniati et al.||2001nRicker & Sarazinn2001> .Subramaniani 



et al. 20 06|) cap able of amplifying magnetic fields by the dynamo 
process ( Subramanian et al. 2006). 

A potential application of small-scale dynamos concerns the 
formation of the first stars and the first galaxies where high- 
resolution simulations have revealed the ubiquity of turbulence in 
the early Universe suggesting that the primordial gas is highly tur- 
bulent ((Abel et al. 2002, O'Shea & Norman 2007, Wise & Abel| 
|2007l[Yoshida et al. 2008 , Greif et al..2008, Turk et al. 20091 ). This 
has strong implications for the exponential amplification of mag- 
netic fields by the small-scale dynamo process during the forma- 
tion of the first stars and in first galaxies ( Schleicher et al.|| 2010| 
|Sur et al.||20T0l [Federrath et al.| |2011b; T urk et al.| |2011). High- 
resolution three-dimensional magnetohydrodynamics (MHD) sim- 
ulations of collapsing magnetized primordial clouds by |Sur et al.] 
([20Tm and'F ederrath et al.|pOTT bl) (hereafter Paper I and II re- 
spectively) show that the turbulence is driven by the gravitational 
collapse on scales of the order of the local Jeans length. This leads 
to an exponential growth of the magnetic field by random stretch- 
ing, folding and twisting of the field lines. In the kinematic stage, 
the dynamo amplification occurs on the eddy-turnover time scale, 
ted = l/v where, / and v are the typical turbulent scale and turbu- 
lent velocity respectively. In a collapsing magnetized system, field 
amplification by the turbulent dynamo should occur on a time scale 
smaller than the free-fall time, tff ^ 1/ \/G p to enable field growth 
faster than the rate at which gravitational compression would am- 
plify the field (Papers I and II). Here p is the total mass density and 
G is the gravitational constant. We note here that the efficiency of 
the dynamo process depends on the Reynolds number and is thus 
related to how well the turbulent motions are resolved in numerical 
simulations. Higher Reynolds number would yield faster field am- 
plification. Indeed, in the limit of infinite magnetic P randtl number 
(i.e., Pr = Rm/Re oo) p] [Schober et al.| ( [20TT| find a depen- 
dence of the growth rate on Re, 



/Re 



(1-0/(1+0 



(1) 



where ^ = 1/3 and / = 1.027 or ^ = 1/2 and / = 0.184 for 
Kolmogorov and Burgers turbulence respectively. 

However, the exorbitant computational costs associated with 
these simulations and the fact that current simulations largely un- 
derestimate the growth rate of the dynamo due to the modest 
Reynolds numbers achievable, restrict following the field evolution 
to only the kinematic phase of the dynamo. An important ques- 
tion concerning the magnetic field growth by the small-scale dy- 
namo process is at what strength does the dynamo generated field 
saturate and what is the structure of those fields. Addressing this 
issue is important to obtain an estimate of the typical saturated 
field strengths to be expected in the protostellar cores. Saturation 
of the field is expected to occur when back reactions either via 
the Lorentz force ( |Subramanian||1999| [Schekochihin et al.|[2004l ) 
or via non-ideal MHD effects such as ambipolar diffusion ( [Pinto | 



|& Galli|[200 8^ become important. In simulations of dynamo ac- 
tion with randomly forced turbulence in a box, saturation of the 
dynamo is achieved when the magnetic energy associated with the 
small-scale field grows to a fraction of the e quipartition value ((Hau 



gen et al.|2 004, Schekochihin et al.|2004[ [Brandenburg & Subra- 



manian 2005( |Federrath et al.| |2011a). However, the exact mean- 
ing of saturation and saturated field strengths in a self-gravitating 
system is far from clear. In this case, the magnetic field is ampli- 
fied by both turbulence and gravitational compression of the field 
lines. Therefore, saturation in a self-gravitating system involves 
also the gravitational energy and hence the ratio of magnetic to 
kinetic energies ^m/^k, does not converge to a constant value 
i.e., Em/E]^ ^ p^l^ (using the fact that the collapse speed v(^p) 
approaches a constant during the collapse ('L arson[[T969 |). Here, 
^ B^V and pv^ V with B being the magnetic field, 

V, p and V are the velocity, density and the volume respectively. 



Thus, we use B ' 



^2/3 



to define dynamo saturation in our case. 



^ Rm and Re are the magnetic and fluid Reynolds numbers respectively. 



The dynamo amplification of magnetic fields is driven by 
turbulent fluid motions. Detailed parameter study of small-scale 
dynamo action in driven turbulence in a box by [Federrath et al.[ 
( [2011a[ ) show that the growth rate and the saturation level of the 
dynamo are sensitive to the Mach number and the turbulence in- 
jection mechanism. It is therefore crucial to explore how variations 
in the turbulent properties of a gas cloud influence the small-scale 
dynamo evolution and saturation in self-gravitating systems. In this 
spirit, we investigate the effects of environment on the gravitational 
collapse and magnetic field amplification in this paper. The key 
questions that we intend to address are the following: how does the 
collapse, growth and saturation level of the magnetic field depend 
on the initial strength and injection scale of the turbulent velocity? 
What information concerning dynamo saturation can be obtained 
from the spectra of magnetic fields? We also seek to understand the 
collapse dynamics and the field amplification when the imposed 
uniform rotation of the cloud is varied. We address these questions 
using the initial conditions described in Papers I and II focussing 
on the collapse of a primordial gas cloud. 

An important issue when studying such systems is the choice 
of the initial field strength. The field strengths obtained from either 
cosmological processes like Inflation or phase transition mecha- 
nisms or from astrophysical mechanisms are weak and have large 
uncertainties (see [Grasso & Rubinstein|200T] for a review). Recent 
FERMI observations of TeV-blazars by Nero nov & Vovk[ poTO] ) 
and Ta vecchio et al.[ |2010) have however reported a lower bound 
of about (10~^^ - 10~^^) Gauss for the primordial field. Numer- 
ical simulations starting with such weak values for the initial field 
strength render it almost impossible to probe the saturation level 
of the dynamo generated fields due to the ^ Re^^^ scaling of the 
magnetic field growth rate (therefore requiring higher numerical 
resolution) and the high computational costs associated. Therefore, 
it is reasonable to start with a stronger seed field which allows us 
to probe the field amplification in both the kinematic and the non- 
linear stage of the collapse. In all the simulations presented here, 
we therefore start with an initial field strength of ^ 1 pG with a 
ratio Em/Ey: ^ 10~^ well below equipartition. In this sense, our 
simulations are to be viewed as controlled numerical experiments 
focussing on exploring the influence of certain parameters (e.g., the 
initial strength and injection scale of turbulence) on the collapse 
and magnetic field amplification. 

The paper is organized as follows. The numerical setup, initial 
conditions and the analysis methods of our simulations are outlined 
in Section 2. We present the results obtained from each of the three 
different parameter regimes i.e., varying the initial strength and in- 
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jection scale of the turbulence, and the amount of initial rotation 
and in Section 3. Finally, we summarize and discuss the implica- 
tions of our results in Section 4. 



2 METHOD 

To study the complex system involving self-gravity, turbulence 
and magnetic fields, we resort to high-resolution three-dimensional 
MHD simulations of a magnetised collapsing cloud using the 
adaptive-mesh refinement (AMR) technique. 



2.1 Numerical setup and initial conditions 

The basic numerical setup is adopted from the one reported in Pa- 
pers I and IL We focus on the gravitational collapse and magnetic 
field amplification of a dense gas cloud, using a simplified setup, 
where we assume a polytropic equation of state, P oc p^, which 
relates the pressure P to the density p with an exponent F = 1.1 
in the density range p ^ (10"^° - 10"^^) g cm~^. Note that 
this almost isothermal equation of state is a good representation 
of the thermal behavior of the primordial gas at the densities con- 
sidered here ( [Omukai et al. 12005) [Clark et al.|2011| ). The numerical 
simulations presented here were performed with the AMR code, 
FLASH2.5 (F ryxell et al.||2000l ). We solve the equations of ideal 
MHD, including self-gravity with a refinement criterion guarantee- 
ing that the Jeans length. 



Aj = 



Gp 



1/2 



(2) 



with sound speed Cg, the gravitational constant G and the density 
p is always refined with a user defined number of cells. The appli- 
cability of ideal MHD in our simulations depends on the strength 
of the coupling between the gas and magnetic fields in primordial 
clouds. Work by Maki & Susa (2 004|[2007] ) focussing on detailed 
models of magnetic energy dissipation via Ohmic and ambipolar 
diffusion find that the ionization degree is sufficiently high in the 
primordial clouds to ensure a strong coupling between ions and 
neutrals, thereby maintaining perfect flux-freezing. This is specifi- 
cally true for the primordial clouds which we address here. On the 
other hand, we expect a higher ionization degree and thus a more 
idealized situation, in the presence of additional radiation back- 
grounds in the later Universe. However, non-ideal MHD effects 
may eventually become important at very high densities, as sug- 
gested by simulations of contemporary star formation (e.g., [Hen] 
[nebelle & Teyssier 2008,puffin & Pudritz 2009 j . These effects are 
not included in the present calculations, but should be the subject 
of future studies. It is to be noted that since we use ideal MHD, 
the magnetic Prandtl number in our simulations is ^ 0(1). From 
earlier studies, we recall that the dynamo amplification requires a 
threshold resolution of about 30 grid cells per Jeans length (see Pa- 
pers I and II). Simulations performed with a resolution below 30 
grid cells are unable to resolve the dynamo amplification of mag- 
netic fields. We therefore perform simulations resolving the local 
Jeans length with a minimum of 64 cells to a maximum of 128 cells 
to explore the influence of initial conditions and the saturated field 
strengths on the small-scale dynamo generated field. We use the 
new HLL3R scheme for ideal MHD ( Waag an et al.|2011|), which 
employs a 3-wave approximate MHD Riemann solver ( Bouchut| 
let al.|2007|[Waagan|2009||Bouchut et al.|2010| ). The MHD scheme 



preserves physical states (e.g., positivity of mass density and pres- 
sure) by construction, and is highly efficient and accurate in model- 
ing astrophysical MHD problems involving turbulence and shocks 
( |Waaganetal.|2011| ). 

Similar to the studies reported in Papers I and II, we model the 
gas cloud as an overdense Bonnor-Ebert (BE) sphere ( Ebert 19551 



Bonnor 



1956) with a core density of pbe = 4.68 x 10 ^ g cm 



at a temperature of T = 300 K. Since the aim of this study is to 
investigate the influence of initial conditions on the evolution of 
the magnetic field, we choose our initial conditions in a way that 
allows us to perform a controlled numerical experiment. The initial 
random seed magnetic field of 1 pG with a power-law dependence 
P{k) ^ was constructed in Fourier space from the magnetic 
vector potential which automatically guarantees a divergence-free 
magnetic field. The turbulence is also modeled with an initial ran- 
dom velocity field with the same power-law dependence as for the 
initial magnetic field. The parameters of the suite of different sim- 
ulations we perform is summarized in Table [T] The parameters a 
and /3 quantify the ratio of initial turbulent energy and the rotational 
energy to the magnitude of the gravitational energy, respectively. 
Note that extremely high values of a may prevent the cloud from 
collapsing under its own gravity. Therefore, we choose initial tur- 
bulent velocities in the range ^;rms 0.2 Cs — 4 Cs such that a varies 
from a minimum of 9.3 x 10~^ to a maximum of ^ 3.723. As for 
the initial injection scale of turbulence, we focus on two cases, one 
where the initial turbulence peaks on scales of the order of the ini- 
tial Jeans length of the core (/inj = 0.7 A j) and the other where it 
peaks on scales /inj = 0.17 Aj. The initial uniform rotation of the 
cloud is varied from /3 = 0% to 8%. 



2.2 Analysis in the collapsing frame of reference 

To understand the behavior of the system quantitatively, we need to 
follow its dynamical contraction in an appropriate frame of refer- 
ence. First, we note that the physical time scale becomes progres- 
sively shorter during the collapse. We therefore define a dimension- 
less time coordinate r (see Papers I and II), 



which is normalized in terms of the local free-fall time. 



tff(t) = 



37T 



32 Gprnit) 



1/2 



(3) 



(4) 



where pm{t) is the mean density in the central Jeans volume, 
Vj = 47r(Aj/2)^/3. If not otherwise stated, we obtain all dy- 
namical quantities of interest within this contracting Jeans volume, 
which is centered on the position of the maximum density. This 
approach enables us to study the turbulence and magnetic field am- 
plification in the collapsing frame of reference. We also note that 
for runs with different initial conditions, like varying degree of ini- 
tial rotation, it is more meaningful to compare different simulations 
at the same central density rather than at the same r. This is be- 
cause, runs with different initial conditions are in different phases 
of the collapse at any given r. Wherever possible, we therefore 
show plots of physical quantities as a function of the central mean 
density. With this in mind, we now proceed to discuss the effects of 
different initial conditions on the gravitational collapse and mag- 
netic field amplification. 
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vJllliUldLlUll 

Run 


1? f^Qollltl C\T\ 

[#of cells/ Aj] 


Tni ti 1^1 Tiifhiil put 
velocity [-Urms] 


Tnitifil Tnipption 

scale [/int] 


Tnitial tiirhiilpnfp 

XlllUlCll UUl UUlCllV^C 

= £^turb/^grav 


(3 = -Erot/I-Egravl 


R64M0.2rotO 


64 


0.2 cs 


0.7 A J 


9.30 xlO-3 


0% 


R64M0.4rotO 


64 


0.4 cs 


0.7 Aj 


3.72x10-2 


0% 


R64M1.0rotO 


64 


1.0 Cs 


0.7 A J 


0.232 


0% 


R64M2.0rotO 


64 


2.0 Cs 


0.7 Aj 


0.930 


0% 


R64M4.0rotO 


64 


4.0 Cs 


0.7 A J 


3.723 


0% 


R64M1.0rot01j0.17 


64 


Cs 


0.17 Aj 


0.233 


0% 


R128M1.0rotO 


128 


Cs 


0.7 A J 


0.232 


0% 


R128M1.0rot4 


128 


Cs 


0.7 Aj 


0.232 


4% 


R128M1.0rot8 


128 


Cs 


0.7 A J 


0.232 


8% 


R16M1.0rot4 


16 


Cs 


0.7 Aj 


0.232 


4% 



Table 1. Summary of the simulation runs. These fall in three broad categories : Runs with varying amplitude 
of the initial turbulence (runs R64M0.2rotO to R64M4.0rotO) followed by a run (R64M1.0rot01j0.17) with a 
different injection scale of the initial turbulence. Finally, runs with with different amounts of initial rotation (runs 
R16M1.0rot4 to R128M1.0rot8). 



3 RESULTS 

In this section, we present the results obtained from numerical sim- 
ulations of magnetic fields in a collapsing environment. We report 
on the influence of varying the initial strength and injection scale 
of the turbulence and, the amount of initial rotation on the col- 
lapse dynamics and magnetic field amplification of a magnetized 
gas cloud. We note that there are two main time scales in the prob- 
lem - the free-fall time, tfr ^ 1 / ^JG p and the eddy turnover time 
scale, ted = l/v{l). For a spectrum which we adopt here, 
v{l) ^ /^/^ and therefore, ted ^^^^ on scales smaller than the in- 
jection scale. In paper II we showed that the driving scale of turbu- 
lence generated during gravitational collapse is of the order of the 
local Jeans length. Klessen & Hennebelle| ( |2010] ) show that it is the 
very process of formation of structures in the Universe on all scales 
that drives the turbulence. Given our initial core density, (subsec- 
tion 2.1) and using the fact that turbulence is driven on the local 
Jeans length (Aj ^ 1.5pc at r = 0) in our simulations, we have 
ted ~ 68.157yr which is much smaller than the initial free-fall 
time scale tff = 0.56 Myr. On a scale ^ 0.1 A j, the eddy turnover 
time scale is even smaller which demonstrates that significant dy- 
namo action will occur on smaller scales provided such scales are 
resolved. In addition to the above two time scales, in presence of 
rotation, there is an additional time scale trot ~ where Q is 
the angular velocity. For /3 = 8%, trot ^ 1.38 Myr. 



3.1.1 Evolution of the density and the velocity 

In order to focus on the effects of turbulence, we do not include any 
ordered rotation of the primordial cloud. Figs.[T]and[2]show the plot 
of different physical quantities against the normalized time, r and 
the mean density pm respectively, for five different values of the ini- 
tial rms turbulent velocity, frms = 0.2cs,0.4cs,Cs,2cs and 4cs. 
The a parameter corresponding to these values are shown in Ta- 
ble [T] We first focus on the hydrodynamic aspects of these curves, 
namely the evolution of the mean density (panel c in Fig. [T| and 
the rms velocity (panels d and c in Figs. [T] and [2] respectively). It 
turns out that initializing the collapse with different values of the 
rms turbulent velocity has an effect on how early or late runaway 
collapse sets in. This is partly due to the fact that the turbulent ki- 
netic energy density prevents the gas from collapsing at the outset 
for initial transonic and supersonic velocities, while for subsonic 
velocities, the gas goes into collapse from r ^ 1. This is shown in 
the evolutionary behavior of the mean density (pm). The run with 
an initial ^;rms ^ 4 Cs does not undergo runaway collapse till about 
r ^ 11 — 12 because the initial a parameter (a = 3.72) is high 
enough to prevent the runaway collapse till the time, the turbulence 
has decayed. The early stage of the evolution is marked by frequent 
changes in the rms velocity (panel (c) in Fig.|2]) in the density range 
(~10~^° — 10~^^)g cm~^. The evolution of the rms velocity for 
runs with initial Vrms = 1 Cs , 2 Cs is marked by a distinct decay 
phase till r ^ 3, after which both curves start to increase as the 
collapse regenerates the turbulence. 



3.1 Effect of initial turbulence 

In this section, we explore the influence of varying the initial rms 
turbulent velocity of the cloud on the collapse and magnetic field 
amplification. The ubiquity of turbulence in the cloud core is cen- 
tral to the idea of magnetic amplification by small-scale dynamo 
action. While subsonic turbulence correspond to typical velocities 
found in the first star-forming minihalos, supersonic velocities may 
correspond to more massive systems like the first galaxies l [Wise| 
let al.|2008||Greif et al.|2008l ). 



3.1.2 Magnetic field evolution and signatures of dynamo 
saturation 

The total magnetic field is amplified by more than 6 orders of 
magnitude, reaching peak values of ^ 1 G at densities of ^ 
10~^^g cm~^ evident from panel (a) in Fig. [2] The evolution 
of the dynamo generated magnetic field is shown in panel (b) of 
Figs. [1] and [2] Except for the run with ^;rms = 4cs, the exponen- 
tial amplification of the magnetic field by the small-scale dynamo 
continues till r ^ 10, illustrating the kinematic phase of the dy- 
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namo. This phase corresponds to densities up to ^ 10~ g cm~ 
for runs with initial ^;rms = (0.2 — 2.0) Cg. The exponential phase 
lasts till pm ^ 10~^^g cm~^ for v^ms = 0.1 Cg. The curves attain 
a peak value of Brms/ Pm ^ in the range ^ (2 — 6) x 10^. The be- 
havior of the run with ^;rms = 4 Cg is different. Here the small-scale 
dynamo field undergoes rapid fluctuations when the core density 
is in the range (10~^° — 10~^^)g cm~^ attaining a peak value 
of 2 X 10^ with respect to the initial value. This is an order of 
magnitude higher than the other runs. It is only after attaining core 
densities of^ 10~^^g cm~^ by when the rearrangement is com- 
plete, that it begins to resemble the other runs. In Fig. [3] we show 
the estimated growth rate of Brms/pm^ computed in the density 
range pm (10~^^ — 10~^^) g cm~^ for runs with different initial 
Mach numbers. This approach ensures that the estimated growth 
rates correspond to the same stage of collapse for all the runs. The 
dynamo growth rate at first increases for subsonic Mach numbers 
and then decreases in the supersonic regime. Table |2] presents a 
summary of the estimated growth rates. Simulations of small-scale 
dynamo action by |Federrath et al.| p011a| ) using external forcing 
also find a similar trend. They find a general increase of the growth 
rate with the Mach number, but a clear drop of the growth rate 
at the transition from subsonic to supersonic turbulence. They at- 
tribute this drop in the growth rate to the formation of shocks at 
transonic speeds, which destroy some of the coherent small-scale 
magnetic field structures necessary to drive dynamo amplification. 
The exponential phase is also marked by an increase in the ratio of 
the magnetic to kinetic energies, Em/Ey, shown in panels (e) and 
(d) in Figs.[l]and[2] 

In the density range pm (lO"^'^ — 10~^°) g cm~^, panel 
(b) in Fig. [2] shows a clear change in the slope of Brms/pm^- As 
we argued in Section [l] such a change of slope with Brms/ pU^ 
either remaining constant or decreasing is an indication of the fact 
that the dynamo proceeds to the saturation phase. The total field 
amplification in this phase is then only due to gravitational com- 
pression of the field lines. For an initial ^;rms = 4.0 Cg, the satu- 
ration phase begins as early aspm 10~^^g cm~^. Panels (e) 
and (d) in Figs.[T]and[2]respectively also show a gradual transition 
to saturation with values reaching about 0.2 — 0.4. Comparing the 
evolution of Brms/pm^ to that of Em/E]^, we find the former to 
be better suitable as a saturation indicator for small-scale dynamo 
action in self-gravitating systems. While there is a distinct change 
in the slope of Brms/pm^ after pm = 10~^^ g cm~^ (panel b in 
Fig. [2]), Em/E]^ still continues to increase mildly attaining values 
in the range 0.2 — 0.4. However, for the run with Vrms = 4.0 Cg, 
both Brms/ Pm ^ and Em/E\^ traces out the saturation phase equally 
well. A convincing probe of dynamo saturation is to explore the 
time evolution of the spectra of the magnetic field. However, since 
the resolution of all the simulations discussed in this subsection is 
only 64 grid cells per Jeans length, we defer the discussion of the 
magnetic spectra to subsection 3.3. Nevertheless, the results ob- 
tained from our Mach number study demonstrate that our choice of 
the initial field strength allows us to probe two regimes of field am- 
plification - an initial phase where both gravitational compression 
and the small-scale dynamo amplify the field followed by a phase 
of saturation of the dynamo in which only compression drives the 
field amplification. 

3.2 Effect of the injection scale 

In this subsection we probe the effect of different injection scales 
of the initial turbulence on the gravitational collapse and magnetic 



Simulation Run 


Mach Number [-Urms/cg] 


Growth rate 


R64M0.2rotO 


0.2 


0.09 


R64M0.4rotO 


0.4 


0.15 


R64M1.0rotO 


1.0 


0.13 


R64M2.0rotO 


2.0 


0.08 



Table 2. Comparison of the growth rates of Brms/Pm ^ as a function of 
the mean density for simulations with different initial Mach numbers taken 
from figure 3. The growth rates are computed when the different runs are in 
the same stage of the collapse. The table shows that the growth rate of the 
dynamo generated field at first increases with the increase in initial subsonic 
Mach number but then starts decreasing for supersonic Mach numbers. 



field amplification by the small-scale dynamo. We consider two 
simulations with the injection scale peaked on scales Zinj/Aj = 
0.17 and /inj/Aj = 0.7. As before, we do not include any ordered 
rotation for these runs and use a resolution of 64 cells to resolve the 
local Jeans length. Both the simulations start with the same value 
of the transonic velocity dispersion. 

Figure |4] shows the variation of the various physical quanti- 
ties as a function of the mean density, pm- With a smaller injection 
scale, the turbulence decays faster (panel (c) in Fig.|4]) compared to 
the run where /inj = 0.7 A j as more kinetic energy has to be ini- 
tialised on the smaller scales to get the same overall turbulent Mach 
number. This can also be seen by comparing the density snapshots 
in the first column of Fig. |5] where the arrows denote the velocity 
vectors. These correspond to times when the central core density 
in both the runs are 10~^^ g cm~^. The upper plot in this col- 
umn corresponds to a simulation with /inj/Aj = 0.17 while the 
plot on the lower panel corresponds to /inj /Aj = 0.7. The run with 
a smaller injection scale shows less turbulent motions inside the 
Jeans volume compared to the run where the initial injection scale 
is peaked on scales of the order 0.7 Aj. 

The second column in Fig. |5] shows that the total magnetic 
field attains a peak value of 10~^ G for /mj/Aj = 0.17 
and about ^ 10~^ G for the run with /inj/Aj = 0.7 at pm = 
10~^^g cm~^. The faster decay in the initial turbulence for the 
smaller injection scale leads to the decay in Brms/pm^ as evi- 
dent from panel (b) in Fig. |4] Over time, as the turbulence gets 
regenerated during the collapse, the magnetic field and hence the 
magnetic energy gets amplified by the dynamo process. The kine- 
matic phase of the dynamo extends topm 10~^^g cm~^ for 
/inj = 0.17 Aj with a peak value of Brms/ pU^ ^ 2.5 x 10^ 
and up to ^ 10~^^ g cm~^ for /inj = 0.7 Aj with a peak value 
of ^ 7 X 10^. Saturation of the dynamo is once again clearly illus- 
trated by the change in slope of Brms/ pU^ . For /inj = 0.17 Aj, this 
occurs from a density '-^ 10~^° g cm~^ while for /inj = 0.7 Aj, 
saturation starts from pm ~ 10~^^g cm~^. Comparing this with 
the evolution of Em/E\, in panel (d), we find that the ratio of the 
energies still continues to increase up to pm 10~^° g cm~^ for 
/i^j = 0.17 Aj and up to 10~^^ g cm~^ for /inj = 0.7 Aj. The 
magnetic energy saturates at about 0.1 of the equipartition value for 
this case while for a smaller injection scale, the saturation level is 
at ^ 0.3 of the equipartition value. 
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3.3 Effect of initial uniform rotation 

Finally, we explore the effects of rotation on the collapse and mag- 
netic field amplification. We begin by comparing simulations per- 
formed at the same resolution (128 cells per Jeans length) with the 
initial rotation parameter (3 having values of 0%, 4% and 8%. 

3.3. 1 Evolution of the density and the velocity 

In Fig. [6] we show the evolution of various physical quantities as a 
function of r for runs with (3 — 0%, 4% and 8%. The same physi- 
cal quantities are plotted as a function of pm in Fig.|7] As reported 
earlier in Papers I and II, the dynamical evolution of the system 
shows two distinct phases. First, as the initial turbulent velocity de- 
cays, the system exhibits weak oscillatory behavior (up to a time 
r — 4) with the mean density pm evolving similarly for runs with 
0%, 4% and 8% initial uniform rotation. This is evident from pan- 
els (c) and (d) in Figs.[7]and[6]respectively. After this, runaway col- 
lapse sets in with peak densities of the order of ~ 10~^^ g cm~^ 
being attained for /3 = 0% and 4%. Higher value of the initial rota- 
tion (e.g. 13 = 8%) leads to a delayed collapse as the cloud is more 
rotationally supported compared to the other two cases (panel (c) 
in Fig.|6|. After the initial decay, the rms velocity increases as the 
turbulence is regenerated by the collapse. 

3.3.2 Magnetic field evolution 

Comparing the evolution of the rms magnetic field in the three 
runs we find from panel (a) in Fig.jT] the magnetic field increases 
by more than 5 orders of magnitude starting from 1 pG reaching 
strengths of 0.1 G at pm ^ 10~^^g cm~^. The corresponding 
evolution with respect to r shown in panel (a) of Fig.|6]may seem to 
suggest that the field amplification is weaker for the 8% run com- 
pared to the other two. Panel (b) in both Figs. [6] and |7] shows the 
evolution of the magnetic field arising due to turbulent motions, i.e, 
due to small-scale dynamo action. For all the three different runs, 
the field is exponentially amplified in the initial phase of decaying 
turbulence (till about r = 3 — 4). After this, the dynamo continues 
to amplify the field up to r ^ 13 or Pm '^10 g cm but the 
amplification is no longer exponential. The ratio of the magnetic to 
kinetic energies (panels e and d in Figs.[6]and|7]also show similar 
behavior with an initial exponential increase followed by a phase 
of linear growth. 

3.3.3 Dynamo saturation 

Saturation of the dynamo occurs from pm ^ 10~^^ g cm~^ when 
all the three curves of Brms/ pm^ shows a change in slope. From 
panel (b) in Fig. [7] these curves tend to remain almost constant 
(red curve) or decay slightly (black and blue curves). However, the 
Em/Ey^ ratio plotted in panel (d) continues to increase with a weak 
dependence on the mean density eventually saturating at values of 
0.2 — 0.4. This once again highlights the fact that Brms/pm^ is 
a better indicator of dynamo saturation in gravitating systems. A 
more detailed analysis of dynamo saturation can be obtained from a 
Fourier analysis of the spectra of the magnetic field. We recall that 
in Kolmogorov turbulence, the eddy turnover time ted = l/v ^ 
^2/3 jj^^g^ smaller scale eddies amplify the magnetic field faster 
due to their shorter eddy turnover times. Therefore, saturation of 
the magnetic field should first occur on the smaller scales and then 
gradually be attained on the larger scales. To illustrate this phenom- 
ena in a more quantitative fashion, we perform a Fourier analysis of 



the magnetic field spectra in the collapsing frame of reference for 
one of our rotation runs with (3 = 8%. To compute the spectra, we 
extract the AMR data in a cube about three times the size of the lo- 
cal Jeans length. For more details on the data extraction procedure 
and the Fourier analysis we refer the reader to subsection 2.4 of 
Paper II. Figure [8] shows the time evolution of the magnetic energy 
spectrum as a function of the wavenumber, k/kj, i.e., normalized to 
the local Jeans wavenumber. Consistent with our initial conditions, 
the spectrum at r = shows the P{B) oc k~^ power-law scaling. 
In addition, the initial spectrum is peaked at k/kj ^ 1.4 which 
is in good agreement with our initial conditions where the initial 
Jeans length of the core is 1.5 pc and the peak of magnetic power 
spectrum is at 0.8 pc. The time evolution of the magnetic spectrum 
shows two important features. Firstly, the peak of the initial mag- 
netic spectrum quickly shifts to smaller scales starting from the ini- 
tial k/kj = 1.4 and then stays roughly constant at k/kj ^ 3 — 4 
till r = 12. This corresponds to about 43 - 32 grid cells which 
is consistent with the Jeans resolution criterion proposed in Paper 
II. Next, from r = 14, the peak of the spectrum shifts to smaller 
k/kj values, i.e., to larger scales. The shift in the peak of the spec- 
trum from k/kj ^ 3 — 4 to k / kj ^ 1 at late times implies that 
the magnetic field at first saturates on the smaller scales followed 
by saturation on larger scales. Thus, the information obtained from 
the time evolution of the magnetic field spectra is consistent with 
theoretical predictions ( Subramanian 1999, Brandenbu rg & Sub-| 
|ramanian|2005| ). The magnetic energy however continues to grow 
with time because the field continues to be amplified by the gravi- 
tational collapse. 

To further analyse the level of saturation we plot the time 
evolution spectra of the ratio of the magnetic to kinetic ener- 
gies in Fig. [9] Since we extract about three times the volume for 
our Fourier analysis, we have taken a mean density of pm ^ 
10~^^ g cm~^ resulting in an initial value of Em/Ey, ^ 2 x 10~^. 
We also remove all large-scale velocity contributions (i.e., global 
rotation and infall) from the kinetic energy. Since the smaller scale 
eddies amplify the field faster, Em/E]^ first starts to peak at large 
k/kj for r > 0. At r = 12, the magnetic energy attains equipar- 
tition with the kinetic energy on a scale k/kj = 20. Once satura- 
tion has been attained on this scale, the peak of the spectrum now 
starts to shift to smaller k/kj.By r = 17.4, the peak of the spec- 
trum occurs at /c/Zcj = 7—10, i.e., saturation of the dynamo is 
now achieved on larger scales. This gradual shift of the peak of 
Em/E]^ ^ 1 to smaller k/kj demonstrates the saturation of the 
dynamo and a possible development of a large-scale coherent mag- 
netic field if the simulation is evolved further in time. The spectra 
further show that the magnetic field becomes dynamically impor- 
tant to back react on the flow on scales k/kj = 20 by r = 12 
with the magnetic energy attaining equipartition with the kinetic 
energy. By the end of the simulation, equipartition is achieved on 
even larger scales. 

3.3.4 Radial profiles 

Figure [T0| shows the radial profile of the density, the rms magnetic 
field, the polar and azimuthal velocity dispersions, a^^e , crv,^ and 
the ratio of magnetic to kinetic energies (Em/E]^) for runs with 
/3 = 0, 4 and 8% initial rotation at a time when all the three at- 
tain the same core density. Similar to earlier results in Papers I and 
II, the density develops a flat inner core and falls off as r~^'^ due 
to the effective equation of state with F = 1.1 (see [Larson|1969| ). 
Panel (b) shows that the rms magnetic field attains peak values be- 
tween 0.02 — 0.04 G within the Jeans volume. The radial profile of 
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Brms shows SL radial dependence oc (r~^"^ — r~^"^) for the three 
runs. This is significantly steeper than the expectation from pure 
flux freezing where B{r) oc r~^/^. The velocity dispersion shown 
in panels (c) and (d) increases in the envelope and drops inside the 
Jeans volume. This is most likely due to the back reaction of the 
strong magnetic fields generated inside the Jeans volume. The ra- 
tio of the magnetic to kinetic energies plotted in panel (e) attains 
values in the range ^0.3 — 0.8 inside the Jeans volume. 

3.3.5 Morphological features 

What effect does rotation have on the morphology of the cloud 
and does it lead to any change in the morphology of the small- 
scale dynamo generated field? In Fig.[TT]we show two-dimensional 
snapshots of the density field in two planes: x — y and x — z for 
the run with /3 = 8% at a time when the central core density is 
^ 10~^^g cm~^. The upper panel shows the zoomed-in slices, 
while the lower panel shows the zoomed-out version of the density 
field. Turbulent motions are dominant inside the Jeans volume as 
can be seen from the plotted velocity vectors in the zoomed-in slice 
in the X — z plane. Rotational motion of the gas cloud is clearly seen 
from the X — 2/ slice in the lower panel. In conformity with known 
results, uniform rotation leads to a flattening of the gas cloud shown 
in the x — z slice illustration in the lower panel. 

Quite interestingly, the morphology of the magnetic field in 
the saturated state does not seem to depend significantly on the 
amount of initial rotation (see Fig. [12]). The magnetic field is ran- 
domly oriented within the local Jeans volume. This is reminis- 
cent of the typical magnetic field structures seen in simulations of 
small-scale dynamo action where turbulence is driven artificially 
with random forcing (Brandenburg & Subramanian 2005^ |Feder-| 
[rath et al.|2011a| ). In our case, the turbulence is driven completely 
by the gravitational collapse ( [Klessen & Hennebelle|2010| Papers 
I and II). In Fig. [12] we plot the slices of the total magnetic field 
in X — ^ and x — z planes for two different cases: the first row 
corresponds to no net rotation and the second row corresponds to 
8% initial rotation. In both cases, the magnetic field is randomly 
oriented, attaining peak values of 0.1 /xG within the central Jeans 
volume. In Fig. [13] we show the radial profiles of the toroidal veloc- 
ity (left panel) and the ratio of the Alfven velocity normalized 
to the toroidal velocity (right panel) at different times. The toroidal 
velocity increases initially as we move from the outer to inner radii 
and then drops inside the local Jeans radius. This is because, in- 
side the Jeans volume, a protostellar core of uniform density (see 
panel (a) in Fig. [To] where the density forms a flat inner core inside 
the Jeans radius) has formed which rotates like a solid body. The 
reason why uniform rotation does not introduce any morphological 
change in the magnetic field can be explained with the help of the 
radial profile of the ratio of the Alfven to the rotational velocity 
shown in the above figure. At any given time, the increase in the ra- 
tio VA I V(f) as one proceeds from the outer to inner regions implies 
that the magnetic field is amplified on a much shorter timescale 
than the orbital time of the cloud. The dynamical effect of rotation 
on the magnetic field morphology will become observable at much 
later times, when the orbital timescale and the timescale on which 
the magnetic field is amplified become comparable. 

3. 3. 6 Resolution effects 

To make a direct comparison of how the collapse and field am- 
plification are affected by rotation, we compare two simulations 



with the same value of the /3 parameter, but at two different reso- 
lutions, one at 16 cells and the other at a resolution of 128 cells. 
The utility of this approach lies in the fact that when the local Jeans 
length is resolved by only 16 cells, the small-scale dynamo is not 
excited (as turbulent motions are under-resolved) and the magnetic 
field is only amplified by gravitational compression in a rotating gas 
cloud. The comparison is shown in Fig. [14] between a 16 cell and 
a 128 cell run both having (3 — 8% initial uniform rotation. For the 
16 cells run, the field amplification arises purely due gravitational 
compression of the field lines. This leads to an order of magni- 
tude difference in the evolution of the rms magnetic field (panel (a) 
in Fig. Q4]). However, as explained earlier, the growth rate of the 
small-scale dynamo depends on the resolution. Thus, an order of 
magnitude difference in this case should be taken as a lower limit. 
The plot in panel (b) confirms the absence of any small-scale dy- 
namo for the 16 cell run as the Brms/pm^ curve at first decays 
and then stays roughly constant while for the 128 cell run, there is 
an initial increase in Brms/pm^ followed by saturation at a value 
Brms/pm^ ^ 6 X 10^. Also, from panel (d), we find that while 
the ratio of Em/E]^ increases rapidly to values ^ 10~^ for the 128 
cell run, there is no corresponding increase for the 16 cell run. We 
refer the reader to Papers I and II for a more detailed analysis of 
resolution effects. 



4 SUMMARY AND CONCLUSIONS 

In this paper, we have presented a detailed study of the influence of 
initial conditions on the gravitational collapse and magnetic field 
amplification of a dense gas cloud. We chose initial and environ- 
mental conditions that are reminiscent of the conditions in primor- 
dial mini-halo s that lead to the formation of the first stars in the 
Universe. We purposely chose the initial strength of the magnetic 
field to be ^ 1 /xG with an Em/Ey, ~ 10~^ to capture the satura- 
tion of the small-scale dynamo. In this respect, this study goes be- 
yond our earlier studies (Papers I and II) where we only considered 
the kinematic regime of the dynamo growth. Since dynamo ampli- 
fication of the magnetic field occurs due to turbulent fluid motions, 
we varied the strength and injection scale of the initial turbulence in 
this study. In addition, we have also studied the influence of initial 
uniform rotation on the collapse and magnetic field amplification. 
The general behavior of the magnetic field amplification presented 
here reveals two distinct phases: first, the total magnetic field is am- 
plified by both the gravitational compression and the small-scale 
dynamo. Later, the small-scale dynamo saturates after which the 
field amplification is driven only by compression. Our main results 
obtained from the systematic study are summarised as follows: 

• We started by exploring the effect of varying the initial 
strength of the turbulent velocity. The parameter range we explored 
concerns an initial Vrms/cs = 0.2,0.4,1.0,2.0 and 4.0. The 
total magnetic field in these simulations is amplified by more 
than 6 orders of magnitude, reaching peak values of ^ 1 G at 
densities of '-^ 10~^^g cm~^. Initially, the small-scale dynamo 
provides additional field amplification to the regular amplification 
occurring due to gravitational compression. Later, we find that the 
small-scale dynamo saturates. In almost all the cases, except when 
the initial ^;rms = 4cs, the amplification of the magnetic field by 
the small-scale dynamo continues till ^ r = 10. This corresponds 
to the kinematic phase of the dynamo. In this phase, Brms/pm ^ at- 
tains a peak value in the range ^ (2 — 6) x 10^ (with respect to the 
initial value) in the density range pm ^ (10~^^ — 10~^°)g cm~^. 
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The growth rate of Brms/pm ^ at first increases with the increase in 
initial subsonic Mach numbers and then decreases in the supersonic 
regime. [Federrath et al.| ( |2011a| , also report a similar behavior 
with Mach number in forced turbulence-in-a-box simulations 
of small-scale dynamo action. Beyond pm ^ 10~^^gcm~^, 
Brms/pm^ shows a distinct change in slope indicating that the 
dynamo reaches saturation. The total field amplification in the 
saturated phase is now only due to gravitational compression. The 
magnetic energy still shows a mild increase reaching values of 
0.2 — 0.4 of the equipartition strength for the parameter range we 
explore. Thus, Brms/pm^ is a more suitable indicator of dynamo 
saturation in self-gravitating systems. 

• Next, we explored the effect of varying the initial injec- 
tion scale of turbulence. We compared two simulations with 
hnj = 0.17 Aj and 0.7 Aj. With a smaller injection scale, the 
turbulence decays faster as more kinetic energy has to be initial- 
ized on the smaller scales to get the same overall turbulent Mach 
number. The smaller scales are therefore dissipated more quickly 
resulting in the early collapse by r ^ 1 compared to the other run 
where runaway collapse sets in at r ^ 4. The kinematic phase of 
the dynamo extends to pm 10~^^g cm~^ for /inj = 0.17 Aj 
and up to ^ 10~^^g cm~^ for /mj = 0.7 Aj. Saturation of 
the dynamo is once again clearly illustrated by the change in 
slope of Brms/pm^ comparcd to the evolution of Em/E]^- For 
/inj = 0.17 Aj, this occurs from a density ^ 10~^° g cm~^ while 
for /inj = 0.7 Aj, saturation starts from pm ^ 10~^^g cm~^. 
By the end of the simulation, the magnetic energy attains values 
in the range 0. 1 — 0.3 of the equipartition values for these two runs. 

• Finally, we explored the effect of uniform initial rotation. We 
considered three simulations where the rotation parameter /3 has 
values of 0%,4% and 8%. The dynamical evolution of the sys- 
tem shows two distinct phases: an initial turbulent decay phase 
followed by a runaway collapse phase where the turbulence gets 
regenerated by the gravitational collapse. The dynamo generated 
magnetic field Brms/pm^ saturates for all the three runs when the 
central core density pm 10~^^ g cm~^. The saturation behavior 
is clearly revealed in the time evolution of magnetic field spectra 
(Fig. |8]) where the peak of the spectrum gradually shifts to larger 
scales (i.e., smaller k/kj values) at late times. Consistent with ear- 
lier theoretical predictions, magnetic energy attains equipartition 
with the kinetic energy first on smaller scales and then gradually 
on larger scales. This is shown in Fig.|9] Our resolution study fur- 
ther shows that the dynamo amplification leads to an order of mag- 
nitude difference in the magnetic field amplification. Since the dy- 
namo amplification is resolution dependent (higher resolution leads 
to stronger field amplification), this difference should be taken as a 
lower bound. 

In summary, using idealized simulations, our systematic study 
has probed for the first time, the effect of different initial conditions 
on the dynamo amplification of the magnetic field and its saturation 
for the specific case of primordial cloud collapse. Our choice of the 
initial field strength allowed us to probe both the kinematic and the 
saturation regime of the small-scale dynamo. Detailed analysis of 
the magnetic field spectra shows that saturation is achieved initially 
on smaller scales and then later on, the magnetic field saturates on 
larger scales. In all our simulations presented here, the magnetic 
energy continues to grow at the expense of the kinetic energy even- 
tually attaining values in the range, Em/ 0.1—0.4. This range 
of values are roughly in agreement with previous analytical ( |Sub-| 



|ramanian|1999| ) and numerical work on small-scale dynamo action 
in forced turbulence (Haugen et al. 2 004^ [Brandenburg & Subra-| 
[manian 2005, Federrath et al. 201 la|. For the case of a primor- 
dial cloud collapse which we explore here, the small-scale dynamo 
comes across as an efficient process to generate strong seed fields 
on scales comparable to the scale of turbulence (in our case, the 
local Jeans length). These strong seed fields could potentially lead 
to the generation of large-scale coherent magnetic fields on scales 
much larger than the scale of turbulence (in our case on scales 
k/kj ^ 1) via large-scale dynamo mechanisms ( [Ruzmaikin et al. 
T988", Hana sz et al.|20 09','Jeda mzik & Sigl| 201 1 ). The imphcations 
of such coherent magnetic fields has been previously explored in 
the works of [Pudritz & Silk| ( [T9891 ) and |Tan & Blackman] ( [20041 ). 
Starting from seed magnetic fields such as the ones generated by 
a small-scale dynamo. Tan & Blackm an (2004 ) show that coherent 
magnetic fields can be produced in situ in primordial star forming 
disks. 

Numerical simulations with an initially coherent magnetic 
field have previously shown the formation of jets and bipolar out- 
flows (e.g., Banerjee & Pudritz 2006, Machida et al. 2006, Baner^ 
jee & Pudritz 2007, Machida et al. 2008 ). Recent work by jPietar- 
,ila Graham et al., ( ,201 Oj find evidence of small-scale dynamo ac- 
tion in the solar surface. [Silk & Langer|p006"] ) infer that magnetic 
fields can be amplified by the magneto-rotational instability (MRI) 
in the disk leading to coherent magnetic fields. The presence of a 
small-scale dynamo can potentially further amplify the seed fields 
in the disk. It remains to be seen what effect, if such large-scale 
magnetic fields generated during the primordial collapse, have for 
example - on the formation process of the first stars. Recent hy- 
drodynamical calculations of first star formation ( [Greif et al.|2011[ 
jClark et aL]|2011| ) have shown that the disks that formed around 
the first young stars were unstable to gravitational fragmentation, 
possibly leading to the formation of small binary and higher-order 
systems. Whether or not, magnetic fields generated and amplified 
by dynamo processes can influence this scenario is an open ques- 
tion. Magnetic fields may also have a bearing on the rotation speed 
of the first stars ( [Stacy et al.|201 1i Moreover, magnetic fields gen- 
erated by dynamo processes in the primordial universe and ejected 
in outflows will have implications for the formation of the second 
generation of stars and the first galaxies. More detailed investiga- 
tions including other relevant physics and additional processes like 
primordial chemistry and cooling, non-ideal MHD and radiative 
feedback effects are needed in the near future to address these pos- 
sibilities. 



ACKNOWLEDGEMENTS 

S.S. thanks the German Science Foundation (DFG) for finan- 
cial support via the priority program 1177 Witnesses of Cos- 
mic History: Formation and Evolution of Black Holes, Galaxies 
and their Environment (grant KL 1358/10). S.S further thanks 
lUCAA and RRI for support and hospitality where this work 
was completed. C. F. received funding from the Australian Re- 
search Council (grant DPI 10102191) and from the European Re- 
search Council (FP7/2007-2013 Grant Agreement no. 247060). 
D. R. G. S thanks for funding through the SPP 1573 (project 
number SCHL 1964/1-1) and the SFB 963/1 Astrophysical Flow 
Instabilities and Turbulence. C.F., R.B., and R.S.K. acknowl- 
edge subsidies from the Baden-Wiirttemberg-Stiftung (grant P-LS- 
SPII/18) and from the German Bundesministerium fiir Bildung und 
Forschung via the ASTRONET project STAR FORMAT (grant 



© 2011 RAS, MNRAS 000,[I]|2l 



Magnetic field amplification during gravitational collapse 9 



05A09VHA). R.B. acknowledges funding by the Emmy-Noether 
grant (DFG) BA 3706. Supercomputing time at the Forschungszen- 
trum Jiihch via grant hhdl42 and hhd20 and at the Leibniz Rechen- 
zentrum via grant no. pr42ho and pr321o is gratefully acknowl- 
edged. The software used in this work was in part developed by 
the DOE-supported ASC/ Alliance Center for Astrophysical Ther- 
monuclear Flashes at the University of Chicago. 



REFERENCES 

Abel, T., Bryan, G. L., & Norman, M. L. 2002, Science, 295, 93 
Banerjee, R., & Pudritz, R. E. 2006, ApJ, 641, 949 
— . 2007, ApJ, 660, 479 

Beck, R. 2004, Astrophysics & Space Science, 289, 293 
Beck, R., & Hoernes, R 1996, Nature, 379, 47 
Beck, R., Poezd, A. D., Shukurov, A., & Sokoloff, D. D. 1994, 
A&A, 289, 94 

Bernet, M. L., Miniati, F, Lilly, S. J., Kronberg, R R, & 
Dessauges-Zavadsky, M. 2008, Nature, 454, 302 

Boldyrev, S., & Cattaneo, F. 2004, Physical Review Letters, 92, 
144501 

Bonnor, W. B. 1956, MNRAS, 116, 351 

Bouchut, F, Klingenberg, C, & Waagan, K. 2007, Numerische 

Mathematik, 108, 7 
— . 2010, Numerische Mathematik, 115, 647 
Brandenburg, A., & Subramanian, K. 2005, Phys. Rep., 417, 1 
Bryan, G. L., & Norman, M. L. 1998, ApJ, 495, 80 
CarilU, C. L., & Taylor, G. B. 2002, Ann. Rev. Astron. Astrophys., 

40,319 

Clark, P C, Glover, S. C. O., Smith, R. J., Greif, T. H., Klessen, 

R. S., & Bromm, V. 2011, Science, 331, 1040 
Dolag, K., Bartelmann, M., & Lesch, H. 1999, A&A, 348, 351 
Duffin, D. F, & Pudritz, R. E. 2009, ApJ, 706, L46 
Ebert, R. 1955, Zeitschrift fur Astrophysik, 37, 217 
Federrath, C, Chabrier, G., Schober, J., Banerjee, R., Klessen, 

R. S., & Schleicher, D. R. G. 2011a, Physical Review Letters, 

107, 114504 

Federrath, C, Sur, S., Schleicher, D. R. G., Banerjee, R., & 

Klessen, R. S. 2011b, ApJ, 731, 62 
Fryxell, B., Olson, K., Ricker, P., Timmes, F. X., Zingale, M., 

Lamb, D. Q., MacNeice, R, Rosner, R., Truran, J. W., & Tufo, 

H. 2000, ApJS, 131,273 
Govoni, F, & Feretti, L. 2004, International Journal of Modern 

Physics D, 13, 1549 
Grasso, D., & Rubinstein, H. R. 2001, Phys. Rep., 348, 163 
Greif, T., Springel, V., White, S., Glover, S., Clark, P, Smith, R., 

Klessen, R., & Bromm, V. 2011, ArXiv e-prints 
Greif, T. H., Johnson, J. L., Klessen, R. S., & Bromm, V. 2008, 

MNRAS, 387, 1021 
Hanasz, M., Woltanski, D., & Kowalik, K. 2009, ApJL, 706, L155 
Haugen, N. E., Brandenburg, A., & Dobler, W. 2004, Phys. Rev. 

E, 70, 016308 
Hennebelle, P, & Teyssier, R. 2008, A&A, 477, 25 
Jedamzik, K., & Sigl, G. 2011, Phys. Rev. D, 83, 103005 
Kazantsev, A. P. 1968, Soviet Journal of Experimental and Theo- 
retical Physics, 26, 1031 
Klessen, R. S., & Hennebelle, P 2010, A&A, 520, A11+ 
Kulsrud, R. M., Cen, R., Ostriker, J. P, & Ryu, D. 1997, ApJ, 480, 

481 

Larson, R. B. 1969, MNRAS, 145, 271 



Machida, M. N., Inutsuka, S.-i., & Matsumoto, T. 2008, ApJ, 676, 
1088 

Machida, M. N., Omukai, K., Matsumoto, T., & Inutsuka, S. 2006, 

ApJ, 647, LI 
Maki, H., & Susa, H. 2004, ApJ, 609, 467 
— . 2007, PASJ, 59, 787 

Miniati, F, Jones, T. W., Kang, H., & Ryu, D. 2001, ApJ, 562, 
233 

Neronov, A., & Vovk, I. 2010, Science, 328, 73 
Omukai, K., Tsuribe, T., Schneider, R., & Ferrara, A. 2005, ApJ, 
626, 627 

O'Shea, B. W., & Norman, M. L. 2007, ApJ, 654, 66 
Pietarila Graham, J., Cameron, R., & Schiissler, M. 2010, ApJ, 
714, 1606 

Pinto, C, & Gain, D. 2008, A&A, 484, 17 

Pudritz, R. E., & Silk, J. 1989, ApJ, 342, 650 

Ricker, P M., & Sarazin, C. L. 2001, ApJ, 561, 621 

Ruzmaikin, A. A., Sokolov, D. D., & Shukurov, A. M., eds. 1988, 

Astrophysics and Space Science Library, Vol. 133, Magnetic 

fields of galaxies 
Ryu, D., Kang, H., Cho, J., & Das, S. 2008, Science, 320, 909 
Schekochihin, A. A., Boldyrev, S. A., & Kulsrud, R. M. 2002, 

ApJ, 567, 828 

Schekochihin, A. A., Cowley, S. C, Taylor, S. F, Maron, J. L., & 

McWiUiams, J. C. 2004, ApJ, 612, 276 
Schindler, S., & Mueller, E. 1993, Astron. Astrophys., 272, 137 
Schleicher, D. R. G., Banerjee, R., Sur, S., Arshakian, T. G., 
Klessen, R. S., Beck, R., & Spaans, M. 2010, A&A, 522, A115-h 
Schober, J., Schleicher, D., Federrath, C, Klessen, R., & Baner- 
jee, R. 201 1, ArXiv e-prints 
Silk, J., & Langer, M. 2006, MNRAS, 371, 444 
Stacy, A., Bromm, V., & Loeb, A. 2011, MNRAS, 413, 543 
Subramanian, K. 1999, Phys. Rev. Lett., 83, 2957 
Subramanian, K., Shukurov, A., & Haugen, N. E. L. 2006, MN- 
RAS, 366, 1437 

Sur, S., Schleicher, D. R. G., Banerjee, R., Federrath, C, & 

Klessen, R. S. 2010, ApJ, 721, L134 
Tan, J. C, & Blackman, E. G. 2004, ApJ, 603, 401 
Tavecchio, F, Ghisellini, G., Foschini, L., Bonnoli, G., Ghirlanda, 

G., & Coppi, P 2010, MNRAS, 406, L70 
Turk, M. J., Abel, T., & O'Shea, B. 2009, Science, 325, 601 
Turk, M. J., Oishi, J. S., Abel, T., & Bryan, G. 2011, ArXiv e- 

prints 

Vainshtein, S. I. 1982, Zhurnal Eksperimental noi i Teoreticheskoi 
Fiziki, 83, 161 

Waagan, K. 2009, Journal of Computational Physics, 228, 8609 
Waagan, K., Federrath, C, & Klingenberg, C. 2011, Journal of 

Computational Physics, 230, 3331 
Wise, J. H., & Abel, T. 2007, ApJ, 665, 899 
Wise, J. H., Turk, M. J., & Abel, T. 2008, ApJ, 682, 745 
Xu, H., Li, H., Collins, D. C, Li, S., & Norman, M. L. 2009, ApJ, 

698, L14 
— . 201 1, ArXiv e-prints 

Yoshida, N., Omukai, K., & Hernquist, L. 2008, Science, 321, 669 



© 2011 RAS, MNRAS 000,[l]|9] 



10 S. Sur, C. Federrath, D. R. G Schleicher, R. Barter jee & R. Klessen 



10° 



10 



-2 



m 10"^ 
10"^ 



00 



10' 



PQ 



CO 

I 



s 

O 
OX) 



10' 

-10 

-12 

-14 

-16 

-18 

-20 
-22 




^ 10 

10 

10 

10 

10 

^ 10 
10 



e 

o 



Loooo r 

0.1000 r 



1 0.0100 K ' 



0.001 
0.0001 









10 15 20 

T 



Figure 1. Evolution of the dynamical quantities in the central Jeans volume as a function of r for runs with different initial Mach numbers of the turbulence: 
Vrms/cs = 0.2, 0.4, 1.0, 2.0, 4.0. Panel (a) shows the rms magnetic field strength Brms, (b) the evolution of Brms/pm^, showing the turbulent dynamo 
amplification by dividing out the maximum possible amplification due to perfect flux freezing during spherical collapse, (c) the evolution of the mean density 
Pm, (d) the rms velocity ^;rms, and (e) the ratio of magnetic to kinetic energy, Em/E]^- All the simulations correspond to a resolution of 64 cells and do not 
include any initial rotation of the cloud. 
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Figure 3. Growth rate of Brms 
density interval pm ~ 10"-*^^ — 10"-*^^ g cm 
generated magnetic field at first increases with increasing subsonic Mach numbers and then starts to decrease in the supersonic regime, 



0.2 Cs, 0.4 Cs, 1.0 Cs and 2.0 Cg. The growth rates for the different runs are computed in the 
^ when all the different runs are in identical stages of gravitational collapse. The growth rate of the dynamo 
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Figure 4. Same as in Fig. [5] but here we compare two runs with zero initial rotation, initial turbulence with i;rms/cs = 1.0 and Jeans resolution of 64 cells 
for different initial injection scales of turbulence: 0.7 /j (black solid line) and 0.17 /j (red dashed line). When initial turbulence is injected on a smaller scale, 
the turbulence decays faster leading to an initial decay of Brms/pm ^• 



© 2011 RAS, MNRAS 000,[l]|9] 



14 S. Sur, C. Federrath, D. R. G Schleicher, R. Barter jee & R. Klessen 



0.05 



— 0.00 




-0.05 




log(P [g cm ''J) 
■ -16.8 



V// \ \ , 

-'^^^^UM ! ^ ^ ^ ^ ^^^^^ 




-0.05 



0.00 
X [pc] 



0.05 



km/ sec 



0.05 



— 0.00 






log(B [/xG]) 
12.0 



-0.05 0.00 
X [pc] 



0.05 





H 



log(p [g cm ^]) 




0.00 



-0.05 




-0.05 



0.00 
X [pc] 



0.05 



km/ sec 



0.05 



0.00 



-0.05 





log(B [mG]) 
^■3.0 



-0.05 0.00 0.05 
X [pc] 



Figure 5. Two-dimensional slices of the density and the total magnetic field for runs with different initial injection scale at a time when the central core density 
is ~ 10"-*^^ g cm~^. The upper row is for /inj/Aj = 0.17 while the lower row corresponds to /inj/Aj = 0.7. In the density snapshots, the arrows denote 
the velocity vectors, while those in the total magnetic field snapshots denote the direction of the local magnetic field. The run with ^inj/Aj = 0.17 shows less 
turbulent motions inside the Jeans volume compared to the run with ^inj /Aj = 0.7. The magnetic field attains a peak value of ~ 10~^ G for /inj /Aj = 0.17 
and ~ 10~^ G for /inj/Aj = 0.7 for the same central density of pm ~ 10"-*^^ g cm~'^. 
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Figure 6. Evolution of the dynamical quantities in the central Jeans volume as a function of r for runs with different initial rotation. Panel (a) shows the 
rms magnetic field strength Brms, (b) the evolution of Brms/pm^, showing the turbulent dynamo amplification by dividing out the maximum possible 
amplification due to perfect flux freezing during spherical collapse, (c) the evolution of the mean density pm, (d) the rms velocity -Urms, and (e) the ratio of 
magnetic to kinetic energy, Em/E]^- AH the simulations correspond to a resolution of 128 cells per Jeans length. 
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Figure 7. Same as Fig. [6] but now plotted versus the mean density. 



(a) _ 






(d) 



10 



-12 



© 2011 RAS, MNRAS 000,[T]j9| 



Magnetic field amplification during gravitational collapse 17 




Figure 8. Time evolution of the magnetic field spectra for the run with 128 cells per Jeans length and £^rot / 1 ^grav | = 8%. The spectra are plotted against 
the wavenumber normalised to the local Jeans wavenumber. The initial magnetic field spectrum (r = 0) follows a power-law on scales smaller than the 
peak scale, k/kj ~ 1.4 as determined by the initial conditions. After the initial decay of the magnetic energy, the peak of the spectra shifts to k/kj ^3 — 4 
and stays constant in this range up to r = 12. Because of the smaller eddy turnover times, these high k/kj modes (smaller scales) attain saturation first. The 
field in the k/kj < 3 — 4 modes attain saturation beyond r = 12 with the peak of the spectra now shifted to k/kj ?^ 1 by the end of the simulation. In 
addition to the amplification by the dynamo, the magnetic field is also simultaneously amplified by gravitational compression. This leads to an overall rise of 
the spectrum with time. 
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Figure 9. Time evolution of the spectra of the ratio of the magnetic to kinetic energies for the same simulation as in figure 8. As smaller scale eddies amplify 
the magnetic field faster, the peak of the spectrum initially shifts to smaller scales (i.e., larger /c/A;j values). By r ~ 10 — 12, the magnetic energy attains 
equipartition with the kinetic energy on scales fc/fcj ~ 20. Further time evolution shows that the peak now spreads to larger scales (i.e., smaller /c/fcj values) 
with the field attaining equipartition on scales k/k;^ = 7 — 10 at r = 17.4. This plot therefore illustrates the saturation mechanism of the dynamo in a 
self-gravitating system. 
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Figure 10. Radial profiles of the density in panel (a), the rms magnetic field in panel (b), components of the velocity dispersion, a^^Q and crv,0 in panels (c) 
and (d) respectively and the ratio of magnetic to kinetic energies in panel (e) at a time when runs with different initial rotation reached the same core density 
as shown in panel (a). Dashed vertical lines depict the local Jeans radius. Error bars are shown as solid vertical lines (panels c and d) inside the Jeans radius. 
The data are obtained from simulations with Aj resolved by 128 cells. 



© 2011 RAS, MNRAS 000,[l]|9] 



20 S. Sur, C. Federrath, D. R. G Schleicher, R. Banerjee & R. Klessen 



T = 14 



0.005 



0.000 



-0.005 



-0.010 



-0.015 



log(P [g cm ^]) 





-0.020 



14 



^^-^^ / / y / f 
.__^_.^/ / / / / / 

/ / f / r t f 
^ / / r f t t t 
^^/l f f f t 
/^^^^^; f t t \ 
/•^y^^t f f 
/ / 




-0.010-0.0050.000 0.005 0.010 0.015 
X [pc] 



km/ sec 



0.10 



0.05 



0.00 I 



log(p [g cm ^]) 
■ -15. 



. N N V 
V 

^ V 
V N \ \ \ 



-0.05 



-0.10 



i / / / / ^ ^ 

I i / / / / 

I 1 / / / ^ 

\ I / I J / ^ ^ 

\ \ I 1 / / 

\ \ I i j ^ ^ " 

^^^^^^// 1 \ \ \ \ f 
^ ^^^----^ / / 1 \ \ \ ^ 
' __-^-^^/ / 1 \ \ \ 

^ / / / r t t 
/ / / / r t 1 

^ ^ ^ ^ y y 



^ / f f r f 
^ y y / f r r t 



/ t 
/ / 
LJL 



-0.10 



-0.05 



0.00 
X [pc] 



0.05 0.10 




log(P [g cm 
15.0 



I 



-0.010-0.0050.000 0.005 0.010 0.015 
X [pc] 



5 km/sec 



14 



0.10 



0.05 



0.00 



-0.05 



-0.10 




\ \ U I i i i / / / / / ^ ^ 
\\\\UI////// 
\ \ \ \ \ U i / ///-^ 

III f U \ \ V ^ V, 

^ y// f f lit n \ \ \ \ \ 

'/ / / / f ftU \ \ \ \\\ 

:/ / f tU \ \ \ \ \\ \ 
f f f t \ \ \ \ \ \ \\ \ 

. . f ! ! f \ \ \ W \ \ \\ \ 
r.f^t t t /if.t.t.t. \ \ A.\.\iv 



log(p [g cm ^]) 
15 



^ V V V V ^ 

\ N. \ V 



i 



-0.10 -0.05 



km/ sec 



0.00 
X [pc] 



0.05 0.10 



km/ sec 



Figure 11. Morphology of the cloud: Two-dimensional slices of the density for a run with 8% initial uniform rotation through the center of the collapsing 
core at r = 14. The top row corresponds to zoomed-in slices in the x-y plane and x-z plane, while the bottom row corresponds to zoomed-out slices in the 
x-y plane and in the x-z plane. The arrows denote the velocity field vectors, while the circle denotes the local Jeans volume. The zoomed out x-z slice on the 
bottom row shows the gradual flattening of the collapsed BE sphere due to the initial rotation. 
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Figure 12. Magnetic field structures: Two-dimensional slices of the total magnetic field through the center of the collapsing sore at r = 14. The upper row 
shows the x-y and the x-z plane snapshots for a run with no initial uniform rotation. The lower row shows the x-y and the x-z plane snapshots for a run with 
8% uniform rotation. As before, the circles depict the local Jeans volume. The arrows denote the direction of the local magnetic field. 




Figure 13. Radial profiles of the toroidal velocity V(f) (left panel) and the Alfven velocity vp^ normalised to (right panel) at different times for a run with 
8% initial rotation. The vertical lines indicate the local Jeans radius at different times. Runs plotted here are from simulations where Aj is resolved with 128 
cells. 
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Figure 14. Comparison between two runs with initial rotation of £^rot / 1 ^grav | = 8%, where the Jeans length was resolved by 128 cells (solid line) and 16 
cells (dashed line), respectively. In comparison to our highest resolution simulation (128 cells), this leads to an order of magnitude difference in the magnetic 
field amplification. No dynamo amplification occurs in the 16 cell run as is evident from panel b (see also Papers I and II for a more detailed study of the 
resolution effects). 
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